arXiv: 1509.04166vl [astro-ph.EP] 14 Sep 2015 


Radar scattering of linear dunes and mega-yardangs: Application to Titan 


Philippe Paillou®'’*, Benoit Seignovert*^, Jani Radebaugh'^, Stephen Wald 


Universite de Bordeaux, UMR 5804-LAB, 33270 Floirac, France 
^Universite de Reims, UMR 7331-GSMA, 51687 Reims, France 
‘^Brigham Young University, Department of Geological Sciences, Provo, UT 84602, USA 
‘^Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA 


Abstract 

The Ku-band (13.8 GHz - 2.2 cm) RADAR instrument onboard the Cassini-Huygens spacecraft has revealed the richness 
of the surface of Titan, as numerous seas, lakes, rivers, cryo-volcanic flows and vast dune fields have been discovered. 
Linear dunes are a major geomorphological feature present on Titan, covering up to 17% of its surface, mainly in 
equatorial regions. However, the resolution of the RADAR instrument is not good enough to allow a detailed study 
of the morphology of these features. In addition, other linear wind-related landforms, such as mega-yardangs (linear 
wind-abraded ridges formed in cohesive rocks), are likely to present a comparable radar signature that could be confused 
with the one of dunes. We conducted a comparative study of the radar radiometry of both linear dunes and mega- 
yardangs, based on representative terrestrial analogues: the linear dunes located in the Great Sand Sea in western Egypt 
and in the Namib Desert in Namibia, and the mega-yardangs observed in the Lut Desert in eastern Iran and in the 
Borkou Desert in northern Chad. We analysed the radar scattering of both terrestrial linear dunes and mega-yardangs, 
using high-resolution radar images acquired by the X-band (9.6 GHz - 3.1 cm) sensor of the TerraSAR-X satellite. 
Variations seen in the radar response of dunes are the result of a contrast between the dune and interdune scattering, 
while for mega-yardangs these variations are the result of a contrast between ridges and erosion valleys. We tested a 
simple surface scattering model, with parameters derived from the local topography and surface roughness estimates, 
to accurately reproduce the radar signal variations for both landforms. It appears that we can discriminate between 
two types of dunes - bare interdunes as in Egypt and sand-covered interdunes as in Namibia, and between two types 
of mega-yardangs - young yardangs as in Iran and older ones as in Chad. We applied our understanding of the radar 
scattering to the analysis of Cassini RADAR T8 acquisitions over the Belet Sand Sea on Titan, and show that the linear 
dunes encountered there are likely to be of both Egyptian and Namibian type. We also show that the radar-bright linear 
features observed in Cassini RADAR T64 and T83 acquisitions are very likely to be mega-yardangs, possible remnants 
of ancient lake basins at mid-latitude, formed when Titans climate was different. 
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1. Introduction 

Cassini-Huygens mission has been in orbit around Sat¬ 
urn since June 2004. The Ku-band (13.8 GHz - 2.2 
cm) RADAR instrument onboard the Cassini spacecraft is 
a combined radiometer/altimeter/scatterometer/imaging 
radar that has revealed a various and rich surface of Ti¬ 
tan through its optically-opaque atmosphere (Elachi et al. , 
2004). RADAR has allowed the discovery of numerous 
seas, lakes, rivers, cryo-volcanic structures and vast dune 
fields (Elachi et ah, 2005; Lopes et ah, 2007; Lorenz et ah, 
2006; Radebaugh et ah, 2007; Stofan et ah, 2007). Dunes 
are in particular a major landform on the surface of Titan, 
since large dune fields cover more than 10 million km^ in 
equatorial regions. They are typically 1-2 km wide, with 
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1-4 km spacing, up to 150 m-high, and can reach more 
than 100 km in length, being mainly east-west oriented 
and aligned parallel with time-averaged equatorial winds 
(Lorenz et ah, 2006; Radebaugh et ah, 2008; Lorenz and 
Radebaugh, 2009; Le Gall et ah, 2011). Dunes on Titan 
are the linear type as observed on Earth, the latter ones 
being used as analogues to infer Titans dunes morphology 
(Neish et ah, 2010; Radebaugh et ah, 2010; Paillou et al., 
2014). 

The RADAR instrument shows Titans dune fields as 
dark linear features separated by brighter linear features. 
Different qualitative interpretations have been proposed 
for this radar signature: (1) dark lines are sand covered 
interdunes with brighter features caused by specular reflec¬ 
tion over dunes crests, as for the linear dunes observed in 
the Namib Desert (Neish et ah, 2010), or (2) dark lines are 
the smooth dunes with brighter linear features caused by 
rougher interdunes, where bedrock is exposed, as observed 
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in the Great Sand Sea in Egypt (Paillou et ah, 2014). Be¬ 
sides linear dunes, other natural wind-related structures on 
Earth show a comparable morphology: mega-yardangs are 
wind-abraded landforms which develop in Earths drylands 
(Sahara, Middle-East, Central Asia), where winds tend to 
be unimodal in direction (Goudie, 2007). Mega-yardangs 
are composed of alternating linear ridges and valleys cre¬ 
ated by wind erosion and sediment deflation, and are often 
associated with soft deposits of desiccated lake beds. As 
for linear dunes, radar images of mega-yardang structures 
show alternating dark linear features (the erosion valleys) 
and brighter linear features (the ridges): confusing radar 
images of linear dunes and mega-yardangs is then quite 
possible, especially at the 300 m resolution of the RADAR 
instrument. 

Comparative planetology is a powerful approach to help 
understand the geology of remote planetary surfaces. We 
conducted a comparative study of the radar scattering 
of both linear dunes and mega-yardangs, based on rep¬ 
resentative terrestrial analogues. We considered the linear 
dunes located in the Great Sand Sea in western Egypt and 
in the Namib Desert in Namibia, and the mega-yardangs 
located in the Lut Desert in eastern Iran and in the Borkou 
Desert in northern Chad. We analysed and modelled 
the radar signatures of both linear structures using high- 
resolution (18 m) radar images acquired by the X-band 
(9.6 GHz - 3.1 cm) radar of the TerraSAR-X satellite. We 
used a simple surface scattering model, whose parameters 
were derived from the local topography (ASTER Global 
Digital Elevation Map - GDEM, and Shuttle Radar To¬ 
pography Mission - SRTM data) and from surface rough¬ 
ness estimates (Bayesian inversion), which accurately re¬ 
produces the radar signature for both landforms. We were 
able to discriminate between two types of dunes: those 
with bare interdunes in Egypt vs. those with sand-covered 
interdunes in Namibia; and between two types of mega- 
yardangs: young ones in Iran vs. older ones in Chad. 

We applied our understanding of the radar scattering to 
the analysis of radar images obtained during the Cassini 
T8 flyby over the Belet Sand Sea on Titan: we show that 
the linear dunes there are likely to be of both Egyptian and 
Namibian types, contradicting previous studies which pro¬ 
posed single-type dune scenarios. We also show that the 
bright linear structures observed in radar acquisitions dur¬ 
ing Cassini T64 and T83 flybys are very likely to be mega- 
yardangs, possible remnants of lake beds at mid-latitude 
(around 40°N). 

2. Linear dunes and mega-yardangs on Earth 

Linear dunes on Earth are mainly located in arid equa¬ 
torial regions, and constitute more than half of terres¬ 
trial dunes (Rubin and Hesp, 2009). They are formed 
in a context of moderate sand supply, with bimodal wind 
regimes associated with seasonal changes, producing quasi- 
symmetrical and linear dunes in a direction parallel to the 


one of average annual wind (Bristow et ah, 2000). Lin¬ 
ear dunes can reach hundreds of kilometres in length, and 
more than a hundred metres in height, with an interdune 
separation of the order of a couple of kilometres (Besler, 
2008). Such results of sedimentary transport and deposi¬ 
tion processes have been observed on all bodies of the Solar 
System that have an atmosphere: Venus, Earth, Mars and 
Titan (Zimbelman et ah, 2013). 

Also belonging to the wind-related linear features family, 
are yardangs (or mega-yardangs at a regional scale), made 
of parallel ridges separated by narrow valleys. The latter 
are formed by the wind abrasion of cohesive rocks (Goudie, 
2007). Mega-yardangs are also mainly located in hyperarid 
regions on Earth, and form under the condition of a strong 
unimodal wind, that transports sand and gravel eroding 
soft sediments on a time scale of a million years (Gabriel, 
1938; Ehsani and Quiel, 2008). A favourable condition for 
yardangs formation is the presence of soft deposits over a 
harder bedrock, typically observed in areas of ancient lake 
basins. Mega-yardangs were observed on Mars (De Silva 
et ah, 2010) and might also be present on Venus (Greeley, 
1999). 

We selected four terrestrial sites to study radar scatter¬ 
ing of linear dunes and mega-yardangs. For linear dunes, 
we considered the Egyptian side of the Great Sand Sea, a 
large dune field covering 300 x 700 km in eastern Libya and 
western Egypt (Besler, 2008). The typical length, width 
and height of dunes in the Great Sand Sea are in the range 
of values proposed for Titan, the interdune being most of 
the time uncovered and exposing the underlying bedrock. 
Egyptian dunes are mainly composed of pure silicate, with 
a low dielectric constant close to the values reported for 
Titans surface materials (Paillou et ah, 2008). They have 
been used as terrestrial analogues to develop a scattering 
model that was applied to the study of the radar response 
of Titans dunes (Paillou et ah, 2014). The Namib Desert 
in Namibia (Lancaster, 1989; Bristow et ah, 2000) was also 
considered for its linear dune fields. It is an arid area which 
borders the southwestern African coast over hundreds of 
kilometres, with mostly linear sand dunes presenting sand- 
covered interdunes, and also comparable in size and mor¬ 
phology to those on Titan. These dunes have been used 
as a training model for radarclinometry (reconstruction 
of the topography from the radar signature), with results 
applied to the study of Titans dunes (Neish et ah, 2010). 

As regards mega-yardangs, we considered the young 
structures of the Lut Desert in eastern Iran, firstly de¬ 
scribed by Gabriel (1938) and later studied by Ehsani and 
Quiel (2008). The region corresponds to a Pleistocene 
basin with fill deposits (silty clay, gypsiferous sands) and 
is about 150 km long by 50 km wide, with ridges reach¬ 
ing heights of 100 m (Goudie, 2007). We also considered 
the older mega-yardangs of the Borkou Desert in north¬ 
ern Ghad, formed by erosion of sandstone of Paleozoic age 
(Mainguet, 1968). This is a region located between the 
Tibesti and Ennedi mountains, in a corridor of very strong 
unimodal wind oriented to the northeast. Eroded valleys 
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between sandstone ridges are mostly filled with aeolian 
sand and present an average depth around 50 m (McHone 
et ah, 1996). 

We acquired radar scenes of the four study sites using 
the X-band radar of the TerraSAR-X satellite, launched 
in 2007 (Pitz and Miller, 2010). The TerraSAR-X sensor 
presents two interesting characteristics for our study: (1) 
X-band is the shortest wavelength (3.1 cm) among avail¬ 
able orbital radars for Earth observation, close enough to 
the Cassini RADAR (2.2 cm) and thus sensitive to the 
same surface roughness range and (2) its high-resolution 
(18 m) allows the detailed study of the radar signatures 
of landforms. We used the TerraSAR-X ScanSAR mode, 
in HH polarisation, to cover 100 x 150 km areas at an av¬ 
erage resolution of 18 m (multi-look images resampled to 
8.25 m), with a varying incidence angle from the near to 
the far range of the image. Acquisition parameters of the 
four TerraSAR-X scenes are summarized in Table 1. 

The TerraSAR-X radar scenes are shown in Fig. 1, with 
relevant full resolution extracts presented in Fig. 2. Both 
linear dunes of the Great Sand Sea and Namib Desert are 
roughly northsouth oriented and are comparable in size, 
but they show a different radar signature mainly due to the 
interdune properties: the bare interdune in Egypt (Fig. la) 
appears brighter than the sand-covered one in Namibia 
(Fig. lb), because of the rougher surface associated with 
the exposed bedrock, which causes a more diffuse scatter¬ 
ing of the incident radar wave. Full resolution extracts in 
Fig. 2a for Egyptian dunes and Fig. 2b for Namibian dunes 
respectively show some brighter radar returns on the side 
of the dunes facing the radar illumination, due to the com¬ 
bination of a lower local incidence angle (Blom and Elachi, 
1987) and finer dune structures, some of them acting as 
natural corner reflectors (Paillou et ah, 2014). TerraSAR- 
X’s high resolution also allows us to see the darker “back 
side” of the linear dunes, corresponding to a low radar 
return because of a higher local incidence angle. 

Mega-yardangs of the Lut Desert and of the Borkou 
Desert do not present the same geographical orientation 
(Fig. 2c and d), but both appear generally brighter than 
linear dunes: most of the radar scattering is due to the 
wide ridges of the eroded layer, which is rougher than 
sand deposits at X-band. Narrow linear valleys separate 
the ridges and appear darker because they are filled with 
smoother aeolian deposits (some valleys even contain very 
small dunes). The young mega-yardangs of the Lut Desert 
show complex ridges structures alternating with poorly de¬ 
fined valleys (see Fig. 2c), while the older mega-yardangs of 
the Borkou Desert exhibit a better contrast between wider 
and flatter ridges and deeper erosion valleys (see Fig. 2d). 

3. Description and modelling of the radar signa- 
tnre of dunes and yardangs 

We studied the variations of the radar signal across se¬ 
lected dune and yardang structures, at TerraSAR-X full- 
resolution, in order to qualitatively understand the scat- 



Fig. 1. Four TerraSAR-X ScanSAR scenes of linear dunes 
of (a) the Great Sand Sea in Egypt, (b) the Namib Desert 
in Namibia; mega-yardangs of (c) the Lut Desert in Iran 
and (d) the Borkou Desert in Ghad. See Table. 1 for geo¬ 
graphical coordinates, white squares show locations of full 
resolution images presented in Fig. 2. 
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Table 1 

The four ScanSAR scenes acquired using the TerraSAR-X radar, orbit inclination is 97°. 


Site - Country 

Scene centre Lat./Lon. 

Incidence angle range 

Acquisition date 

Scene ID 

Orbit look 

Great Sand Sea Egypt 

24°30'N 26°10'E 

19.5-30.0° 

2013-06-06 

371348705_2 

ascending right 

Namib Desert Namibia 

24°30'S 15°10'E 

24.5 - 34.5° 

2013-12-27 

372338466_1 

ascending right 

Lut Desert Iran 

30°15'N 58°20'E 

29.0 - 38.5° 

2013-06-08 

371348705.4 

ascending right 

Borkou Desert Chad 

19°00'N 19°05'E 

24.5 - 35.0° 

2013-12-08 

372239342.1 

ascending right 



Fig. 2. Full resolution (8.25 m/pixel) extracts of the 
TerraSAR-X scenes: linear dunes of (a) the Great Sand 
Sea in Egypt (incidence angle 23.0°), (b) the Namib Desert 
in Namibia (incidence angle 29.2°); mega-yardangs of (c) 
the Lut Desert in Iran (incidence angle 33.6°) and (d) the 
Borkou Desert in Chad (incidence angle 33.5°). 


tering mechanisms involved, and then to reproduce the 
observed variations with the help of a simple surface scat¬ 
tering model. 

3.1. Radar signatures of linear dunes and mega-yardangs 

Eight representative locations were selected in the Great 
Sand Sea and six in the Namib Desert, in order to derive a 
“typical” radar scattering profile for Egyptian and Namib¬ 
ian dunes. We applied a standard calibration procedure to 
full resolution TerraSAR-X data to compute the backscat- 
tered radar power at each pixel (Schwerdt et ah, 2010). 
We averaged pixels values in the direction parallel to the 
main orientation of the structures, in order to get a better 
representation of typical variations of the radar scattering 
across the studied structures: averaging was performed 
using a median estimator, so that extreme high returns 
due to natural corner reflectors and extreme low returns 
at the noise floor level (around -28 dB for TerraSAR-X) 
did not bias the average value. Fig. 3 shows two typi¬ 
cal radar signatures across Egyptian and Namibian linear 
dunes. Egyptian dunes present three main levels of radar 
scattering (Fig. 3, left): a radar-bright interdune (around 
-15 dB) due to an exposed rough bedrock (Besler, 2008), a 
medium scattering level for the dunes side facing the radar 
illumination (“front side”, around -23 dB), and a low scat¬ 
tering level for the dunes side opposite to the radar illumi¬ 
nation (“back side”, -30 dB and lower). Namibian dunes 
present a different, rather bimodal, radar signature, with 
a less radar-bright interdune due to a smooth sand cover 
(Bristow et al., 2000). Fig. 3 (right) shows a typical radar 
profile across a Namibian dune, where two main scattering 
levels can be observed: a strong radar return on the side of 
the dune facing the radar, can reach -10 dB (due to a lower 
incidence angle combined with smaller-scale structures su¬ 
perposed to the main dunes shape (Lancaster, 1989). A 
lower scattering level, in the -25 to -30 dB range, corre¬ 
sponds to returns from both interdune and back side of the 
dune (the latter likely mixing with the interdune return 
when decreasing the radar image resolution). This shows 
that, while both dune systems in Egypt and Namibia are 
of the same linear type, their radar signature can be quite 
different due to different interdune properties. This ap¬ 
pears quite obvious here from high-resolution radar im¬ 
ages, but at the lower resolution (and higher noise level) 
of the Cassini RADAR instrument, observing these struc¬ 
tures can become more confusing. 

Regarding mega-yardangs, we studied five representa¬ 
tive locations in the Lut Desert and also five in the Borkou 
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Fig. 3. TerraSAR-X radar scattering and topography (SRTM with voids filled using GDEM) profiles across a linear 
dune of the Great Sand Sea in Egypt (left) and of the Namib Desert in Namibia (right). The analysed full resolution 
TerraSAR-X extracts are shown on top of the figure (site #4 for the Egyptian dune, at incidence angle of 23°, and site 
^1 for the Namibian dune, at incidence angle of 29.2°), north is up and radar illumination is from the left. Vertical axis 
of plots is backscattered power (in dB) or altitude (in m) while horizontal axis is the distance across the dunes (in m). 
Location of interdune, front side and back side of the dunes are indicated on each plot. 
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Desert, to derive typical radar profiles for young and older 
mega-yardangs. As before, we averaged the pixel values in 
the direction parallel to the main orientation of the struc¬ 
tures using a median estimator. Fig. 4 presents two typi¬ 
cal radar signatures across Iranian and Chadian yardangs. 
The young mega-yardangs in Iran (Fig. 4, left) show a high 
spatial frequency alternation of the high scattering level 
from the rough crests of the eroded layer (>-10 dB) and the 
medium-scattering level from narrow and shallow valleys 
(around -15 dB). This alternation is cut in places by re¬ 
gions of low radar return (<-25 dB), due to deeper eroded 
valleys filled with smooth aeolian deposits (Gabriel, 1938). 
The older mega-yardangs in Chad present a different, more 
bimodal radar signature (Fig. 4, right), with alternation of 
radar-bright level from wider and flatter crests (around -10 
dB, less radar-bright than the Iranian ones likely because 
a longer period of erosion produced a smoother surface), 
and low-scattering level due to larger erosion valleys (<-25 
dB), which are likely to be filled with aeolian sand deposits. 
As for linear dunes, the two mega-yardangs structures, of 
different age and created in different types of rocks, exhibit 
differences in their radar signature. However, it is appar¬ 
ent that mega-yardangs in general are more radar-bright 
than linear dunes, since we observe about 10 dB difference 
in the average backscattered power. This is likely due to 
the strong radar return of the rougher and wider yardangs 
crests, compared with a lower overall signal returned by 
the smoother sandy surfaces of dunes. 

3.2. A simple model for the radar signature of dunes and 
yardangs 

After presenting a qualitative description of radar scat¬ 
tering profiles across linear dunes and mega-yardangs, we 
tried to reproduce the variations of the observed radar 
backscattered power using simple surface scattering mod¬ 
els, whose parameters were estimated from various sources. 
We considered two surface scattering models, widely used 
by the radar remote sensing community: the Geometric 
Optics Model COM suited to rough surfaces (surface 
roughness being defined with respect to the radar wave¬ 
length, most natural surfaces are rough at X-band) (Fung 
and Eom, 1981), and the Integral Equation Model lEM 
valid for medium-rough to smooth surfaces (Fung et ah, 
1992). The latter model was previously used to reproduce 
the radar scattering of linear dunes of the Great Sand Sea 
in Egypt, using SIR-G/X-SAR scenes (Paillou et ah, 2014). 

Besides the radar wavelength and polarisation, input pa¬ 
rameters for these models are local incidence angle 9 (com¬ 
puted from the radar look angle and local slope at the ob¬ 
served point), surface roughness described by rms-height 
a and correlation length L (assuming a Gaussian auto¬ 
correlation function), and the dielectric constant e of the 
material constituting the surface. The dielectric constant 
of dry sediments and sedimentary rocks in desert regions 
does not vary much and is the 36 range for its real part, 
with an imaginary part close to zero (Ulaby et ah, 1990). 
We fixed the dielectric constant to that of silicate (e = 3.5) 


for all surfaces. Note that small changes in this parame¬ 
ters are not important when computing radar backscat¬ 
tered power, as compared to the effects of variations in 
the incidence angle and surface roughness. The topogra¬ 
phy of the various selected areas was obtained using SRTM 
and GDEM data, GDEM being mainly used to fill in holes 
in the SRTM coverage, especially in dune areas. Taking 
into account the look angle and orbit inclination of the 
TerraSAR-X radar, we computed slope maps in the radar 
range direction, from which we derived the local incidence 
angle 6 at each pixel. The resolution of digital elevation 
models provided by SRTM (90 m) is coarser than the res¬ 
olution of the TerraSAR-X images, so it was not possible 
to reproduce scattering effects due to small-scale features, 
in particular natural corner reflectors. As regards sur¬ 
face roughness parameters, only the ones for the surface of 
Egyptian dunes were estimated by previous studies (Pail¬ 
lou et ah, 2014). We then considered a Bayesian inversion 
approach to estimate the surface roughness of dunes and 
interdunes for both Egyptian and Namibian cases, and of 
crests and valleys for both Iranian and Ghadian yardangs. 
We computed the probability P of similarity between the 
actual radar backscattered power in TerraSAR-X images 
computed radar backscatter power for 
roughness parameters a varying between 0.02 and 0.5 cm 
and L varying between 0.2 and 3.5 cm, assuming a toler¬ 
ance criteria r= 0.03: 


P(cr, L) = - exp 

T 



- ^th L) ^ 


( 1 ) 


The GOM model was used when a > 0.3 cm and the 
lEM model was used when a < 0.3 cm. After explor¬ 
ing the roughness parameter space, we kept the (cr, L) 
couple which corresponds to a maximum probability of 
similarity between the actual TerraSAR-X and the com¬ 
puted backscattered power. Fig. 5 shows an example of 
this Bayesian approach for dunes and interdunes of the 
Great Sand Sea in Egypt. 

Once incidence angle and surface roughness parameters 
were estimated, we computed “theoretical” radar scatter¬ 
ing profiles using GOM and lEM models, and compared 
them to observations for both linear dunes and mega- 
yardangs. Fig. 6 shows examples of comparison between 
computed and observed radar profiles across linear dunes 
in Egypt and Namibia, and mega-yardangs in Iran and 
Ghad. One can see that the computed radar profiles and 
the actual TerraSAR-X ones are quite similar on the aver¬ 
age. The high-frequency variations in TerraSAR-X profiles 
could not be reproduced, due to the coarse resolution of 
the SRTM topography. Nevertheless, simple surface scat¬ 
tering models such as GOM and lEM allow us to explain 
fairly well the variations of the actual radar scattering, as 
previously shown in Paillou et al. (2014) for dunes only. 
These results confirm the qualitative interpretation of the 
radar signature of linear dunes and mega-yardangs pre¬ 
sented in the previous sub-section. 
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Fig. 4. TerraSAR-X radar scattering and topography (SRTM with voids filled using GDEM) profiles across mega- 
yardangs of the Lut Desert in Iran (left) and of the Borkou Desert in Chad (right). The analysed full resolution 
TerraSAR-X extracts are shown on top of the figure (site #1 for Iranian yardangs, at incidence angle of 33.6°, and site 
^5 for Chadian yardangs, at incidence angle of 33.5°), north is up and radar illumination is from the left. Vertical axis 
of plots is backscattered power (in dB) or altitude (in m) while horizontal axis is the distance across the yardangs (in 
m). Location of valleys and crests are indicated on each plot. 
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Fig. 6. Comparison between radar scattering profiles computed using our surface scattering model (bold black line) 
and TerraSAR-X actual data (grey line). From top to bottom: Egyptian dune, Namibian dune, Iranian yardang, and 
Chadian yardang, corresponding to the sites shown in Fig. 3 and Fig. 4. Vertical axis is the backscattered power (in dB) 
and x-axis is a normalised distance across the studied structure. The dash-line represents the TerraSAR-X noise level 
(around -28 dB). 
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4. Application to Titan 
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Fig. 5. Bayesian inversion for the surface roughness of 
an Egyptian dune. The probability P increases from blue 
to red, and the black dot is the location of the maximum 
probability. Top: side of the dune facing the radar illu¬ 
mination (maximum at a= 0.07 cm and L = 1.39 cm). 
Middle: side of the dune opposite to the radar illumina¬ 
tion (maximum at a= 0.09 cm and L = 1.95 cm). Bottom: 
rough interdune (maximum at cr= 0.48 cm and L = 0.63 
cm). The discontinuity in probability at cr= 0.3 cm is due 
to the transition from GOM to lEM model. (Eor interpre¬ 
tation of the references to colour in this figure legend, the 
reader is referred to the web version of this article.) 


Comparative planetology is a powerful approach when 
studying remote bodies, especially when only images of 
their surface are available. This is particularly true for 
Titan, which shows a complex and rich surface where geo¬ 
logical and meteorological processes are active, and where 
fieldwork is still a dream. It is even more difficult to in¬ 
terpret radar images, which are less intuitive than their 
optical counterparts, but we only have such (low resolu¬ 
tion and noisy) radar scenes to “see” the surface of Titan 
through its opaque atmosphere. Studying analogue terres¬ 
trial landforms using a comparable radar system (at least 
in terms of wavelength) can then give us some useful in¬ 
dications and directions to help understand the Cassini 
RADAR images. However, one has to remain cautious, 
since the actual origin, age, morphology, and materials of 
the exotic structures observed on the surface of Titan are 
not well known. With this in mind, the previous qualita¬ 
tive and quantitative study of the radar scattering of ter¬ 
restrial linear dunes and mega-yardangs can help us to in¬ 
vestigate similar linear structures on Titan. We focus here 
on some selected sites on Titan imaged by the RADAR 
instrument, and show how the radar images of terrestrial 
analogues can guide the interpretation. 

Regarding linear dunes, we considered the Cassini 
RADAR acquisition performed during the T8 flyby in Oc¬ 
tober 2005. This covered the more than 2500 km-long 
Belet Sand Sea, a region dominated by linear dunes, and 
stretching from 180°W to 300°W and from 20°N to 20°S 
(Lunine et ah, 2008). Illumination direction is roughly 
perpendicular to the orientation of the main axis of the 
dunes, so that the acquisition geometry is comparable to 
that of TerraSAR-X images presented in the previous sec¬ 
tion. Incidence angle during the T8 flyby varied from 17° 
to 24°, also comparable to the range of incidence angle of 
our TerraSAR-X data. The RADAR instrument resolu¬ 
tion is sufficiently coarse that it can be difficult to deter¬ 
mine if the linear dunes in the Belet Sand Sea have sandy 
interdunes like those in Namibia (Neish et ah, 2010), or 
are separated by exposed bedrock like the Egyptian dunes 
(Paillou et ah, 2014), or if the two kinds can be found 
in different locations across Belet. We studied ten loca¬ 
tions in the T8 flyby, where dune structures are clearly 
imaged, and performed some averaging in the direction 
parallel to the main orientation of the dunes in order to 
extract significant radar scattering profiles. We present 
three interesting and typical radar signatures in Eig. 7. 
Site )(^3 is located in a region where radar-dark dunes end 
on what seems to be a radar-bright plateau, and where 
the extracted radar scattering profile looks similar to those 
observed for Egyptian dunes: a first radar-bright level (- 
4 dB) corresponds to the rough interdune (the exposed 
surface of the plateau), a second intermediate level (-7 to 
-8 dB) might correspond to front sides of the dunes illu¬ 
minated by the radar, and a third lower level (-12 to -13 
dB) could be associated to sides of the dunes opposite to 
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the radar illumination direction. Site ^8 is located in a 
region that was previously studied by Neish et al. (2010) 
who applied radarclinometry techniques to estimate dune 
heights. It shows the same three “Egyptian-like” scatter¬ 
ing levels: a radar-bright level (-1 to -4 dB) that could 
again be associated with bare interdunes, an intermediate 
level (-7 to -8 dB) that could correspond to front sides 
of the dunes, and a lower scattering level (<-12 dB) that 
could be associated with back sides of the dunes (and re¬ 
mains higher than the Cassini RADAR noise level, around 
-18 dB). However, it is also possible that dunes of site 
#8 are of the Namibian type, and that the radar-bright 
peaks are the sides of the dunes facing the direction of the 
radar illumination, the intermediate level being the result 
of interdune scattering, and the lower level being the back 
side of the dunes. Interpreting the radar-bright features 
as interdunes (Egyptian-type) or as sides of dunes facing 
the radar (Namibian-type) will have a strong impact on 
the dunes shape estimation using radarclinometry (Neish 
et ah, 2010). Site #7 exhibits a more clear Namibian- 
like bimodal radar signature. It shows alternation of high- 
scattering return (-6 to -7 dB) and medium-scattering re¬ 
turn (around -12 dB): the signature of site #7 looks very 
similar to that of Namibian dunes, the radar-bright fea¬ 
tures being sides of the dunes facing the radar, and the 
medium-scattering level being related to sand-covered in¬ 
terdunes, mixed with the low return of back sides of the 
dunes. One can also notice that the radar-scattering levels 
for Titans dunes are generally higher than for terrestrial 
ones, possibly due to some additional volume scattering in 
a heterogeneous or porous material, as proposed by Pail- 
lou et al. (2014). While the low resolution (at best 300 
m) combined with the high noise level (around -18 dB) 
of the RADAR instrument does not allow us to extract 
unambiguous radar signatures, both scattering profiles of 
Egyptian and Namibian-type dunes seem to be present in 
the Belet Sand Sea. Eurther studies of Titans dune fields 
should then take into account such possible different in¬ 
terpretations of the dark and bright linear features, and 
evaluate their impact on the understanding of the mor¬ 
phology and history of the structures. 

Concerning possible mega-yardangs on Titan, previous 
work already proposed some candidate structures (Pail- 
lou and Radebaugh, 2013; Radebaugh et ah, 2015). We 
considered two interesting landforms, observed at mid¬ 
latitudes during the Cassini RADAR flybys T64 and T83. 
The T64 flyby was performed in December 2009 (Cassini 
extended mission) and our region of interest is centred on 
41°N/210°W with an incidence angle around 14°. The 
T83 flyby was acquired in May 2012 (Cassini second ex¬ 
tended mission) and our region of interest is centred on 
40°N/197°W with an incidence angle around 12°. Fig. 8 
shows two very bright structures in T64 and T83, pre¬ 
senting alternating dark and bright linear features, which 
are interpreted as possible mega-yardangs. The two land- 
forms are of comparable size, and cover an area of about 
60 X 60 km. Four selected sites were studied in the T64 


structure and three in the T83 one, crossing bright-dark 
alternating linear features (see Fig. 8 presenting radar pro¬ 
files across two representative sites). Again, we averaged 
the radar scattering along the main direction of the linear 
features, in order to derive significant radar scattering pro¬ 
files. While both T64 and T83 candidate mega-yargdangs 
look similar, a closer look shows that the T64 structure is 
more like the old mega-yardangs in Chad, while the T83 
structure is more similar to the Iranian mega-yardangs. 
Fig. 8 presents the radar scattering profile associated with 
site #I in T64 flyby: as for Chadian mega-yardangs, two 
main scattering levels are observed, a very bright level (>0 
dB) corresponding to possible wide yardang crests, and a 
less-bright level (around -5 dB) corresponding to possible 
erosion valleys. Fig. 8 also shows the radar profile as¬ 
sociated with site #1 in T83 flyby: the same two radar 
scattering levels as in T64 are observed there, but with 
a higher spatial frequency alternation of narrower bright 
and less-bright features, more similar to the signature of 
young yardangs observed in the Lut Desert in Iran. We did 
not conduct any morphological measurements on these po¬ 
tential extra-terrestrial yardangs to compare them to the 
terrestrial ones (whose morphology is actually poorly de¬ 
scribed (Goudie, 2007), but we rather focused here on their 
radiometric signature. As for linear dunes, one can notice 
that the structures observed on Titan are much brighter 
than their terrestrial analogues: the brightest crests of Ira¬ 
nian and Chadian yardangs at around -5 dB correspond to 
the low scattering level of the “valleys” of Titans yardangs. 
Again, additional scattering processes, such as volume 
scattering, should be considered on Titan. Nevertheless, 
even if the global radar scattering level on Titans surface 
is higher, the average scattering level difference between 
linear dunes and mega-yardangs remains of the order of 
10 dB, comparable to the difference observed for their ter¬ 
restrial analogues. Finally, it is interesting to note that the 
average radar scattering level of the T64 and T83 struc¬ 
tures is around -2 dB (for an incidence angle around 13°), 
a value comparable to the average scattering level of north 
polar empty lakes basins (Hayes et ah, 2008; Michaelides 
et ah, 2015). If the radar-bright structures in T64 and 
T83 are mega-yardangs, they may be the erosion result of 
unimodal wind direction over soft deposits that could be 
ancient lake beds and/or evaporites (Barnes et al., 2011; 
Cordier et ah, 2013). If confirmed, this would indicate the 
possible existence of ancient lake basins at lower latitudes 
(around 40°N) than polar regions in the past history of 
Titan. 

While comparing terrestrial geological landforms to ones 
observed on Titans surface is not an easy exercise, the 
choice of relevant structures on Earth and the use of a 
high-frequency orbital radar sensor can nevertheless help 
to better interpret the noisy and low resolution images 
provided by the Cassini RADAR instrument. In order 
to illustrate, in a qualitative way, the similarity between 
the radar signature of linear structures on Earth and Ti¬ 
tan, we degraded the four TerraSAR-X scenes that we 
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Fig. 7. Extracts of Cassini RADAR flyby T8 showing linear dunes of sites ^3, #8, and #7 with corresponding 
radar scattering profiles. Site #3 is located 6.5°S/238.8°W with an incidence angle around 21°, site #8 is located 
8.8°S/259.4°W with an incidence angle around 26°, and site #7 is located 8.2°S/248.9°W with an incidence angle 
around 26°. North is down, the “i” arrow indicates direction of the radar illumination. 
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Fig. 8. Extracts of T64 (left) and T83 (right) Cassini RADAR acquisitions showing very bright linear structures, 
interpreted as possible mega-yardangs, together with two studied radar profiles. Resolution is 175 m/pixel, each image 
covers about 100 x 120 km, T64 yardang is located 41.3°N/210.2°W and T83 yardang is located 39.8°N/196.6°W, radar 
incidence angle for both images is around 13°. 
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Fig. 9. The four TerraSAR-X scenes of Fig. 1 with a 
resolution degraded to 300 m and a noise level increased 
to -18 dB. Linear dunes in Egypt (a) and Namibia (b) 
look quite similar to the structures observed in the Belet 
Sand Sea, and mega-yardangs in Iran (c) and Chad (d) are 
recalling the radar-bright structures observed in T64 and 
T83. 

previously used to a resolution and a noise level similar 
to the RADAR instrument. Fig. 9 presents the result¬ 
ing degraded images of linear dunes and mega-yardangs 
on Earth, which are strongly reminiscent of the landforms 
seen on Titans surface. 

5. Conclusion 

High resolution X-band radar images of terrestrial lin¬ 
ear dunes and mega-yardangs can help to better interpret 
the radar signature of similar structures that are likely to 
exist on the surface of Titan. Using a simple surface scat¬ 
tering model combined with topography data, we accu¬ 
rately reproduced the radar backscattering profiles of lin¬ 
ear dunes (Egypt and Namibia) and mega-yardangs (Iran 
and Chad), thus establishing a relationship between the 
radar radiometry and the morphology and nature of im¬ 
aged structures. 

Applying this new understanding of the radar scattering 
process to the interpretation of Cassini RADAR data, we 
have shown that both Egyptian and Namibian-type dunes 
are present in the T8 flyby acquisition of the Belet Sand 
Sea, i.e. interdune regions there can be either sand-covered 
or not. This result should be taken into account for future 
studies of the equatorial dune fields present on Titan. 

We have also shown that two radar-bright structures, 
observed during the T64 and T83 flyby acquisitions, are 


very likely to be the first mega-yardangs observed on Ti¬ 
tan. Such erosional structures could be the remnants of 
past lake basins at mid-latitude, formed when Titans cli¬ 
mate was different, and need for further studies to under¬ 
stand their origin. 

As previously noticed by other authors, both dunes and 
yardangs on Titan present a much stronger radar return 
than their terrestrial counterparts, indicating that some 
additional scattering processes, such as volume scattering 
in an heterogeneous or porous material, occur on Titan. 
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